Optical saturation driven by exciton confinement in molecular-chains: a TDDFT study 
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We have identified excitonic confinement in one-dimensional molecular chains (i.e. polyacetylene 
and H2) as the main driving force for the saturation of the chain polarizability as a function of 
the number of molecular units. This conclusion is based on first principles time-dependent density 
functional theory calculations performed with a new derived exchange-correlation kernel. The failure 
of simple local and semi-local functionals is shown to be related to the lack of memory effects, spatial 
ultranonlocality, and self-interaction corrections. These effects get smaller as the gap of the system 
reduces, in which case such simple approximations do perform better. 

PACS numbers: 71.35.-y, 31.15.oe, 31.15.ap 



The electronic quantum confinement occurring in low- 
dimensional systems is often responsible for the pecu- 
liar spectroscopic properties exhibited by molecules and 
nanostructures. Exciton confinement explains, for exam- 
ple, the size-dependent color of semiconducting quantum 
dots used as fluorescent markers in biology Highly 
localized excitons (solitons) play also a fundamental role 
in describing the process of vibrational energy transfer 
in complex proteins Q. In this context, one-dimensional 
polymers and molecular chains constitute an excellent 
playground to analyze the interplay between correlation 
effects and quantum confinement. For example, in non 
conducting polymers the longitudinal linear polarizabil- 
ity per monomer unit a(N)/N tends to a constant in the 
large-N limit [3| . This optical saturation stems from the 
polarization of the electrons along the chain that tends 
to counteract the external field and manifest itself in the 
dependence of the optical absorption spectra on the poly- 
mer size. The present Letter aims to provide a consistent 
description of this saturation within density-functional 
based schemes. 

The commonly used local (LDA) and semilocal (GGA) 
approximations to Density-Functional-Theory (DFT) 
and Time-Dependent DFT (TDDFT) that successfully 
describe the electronic properties of many different phys- 
ical systems Q , fail dramatically for the case of semicon- 
ducting one-dimensional molecular chains. ALDA do not 
describe the main features of the absorption spectra and 
strongly overestimates the polarizability with respect to 
quantum chemical calculations 0, @, 0] . The reason for 
this poor performance of ALDA is still not settled, and 
it has been traced back to the need of long-range (LR) 
terms (ultra non-locality) in the exchange-correlation 
functionals \f§ or to the lack of self-interaction correc- 
tion (SIC) [a EH. Exact-exchange (EXX) or current- 



density functional (CDFT) approaches capture some of 
the effects, but not all: while EXX works fairly well in 
the case of the finite H2 chain (reproducing the Hartree- 
Fock results 01) it fails in reproducing the absorption 
spectrum of the infinite trans-polyacetylene chain 11 1. 



On the contrary CDFT provides quite good results for fi- 
nite 7r-conjugated polymers [r^ll3j|. while it breaks down 
for the H2 chain. 

In TDDFT all the effects beyond the non interacting 
particles approximation are embodied in the exchange- 
correlation (xc) kernel f xc = 6v xc /6n, with v xc the xc- 
potential and n the exact ground state electronic den- 
sity. The recent developments in merging many-bod^ 
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perturbation-theory (MBPT) [1J| and TDDFT 
open the path to unravel the physical origin of the re- 
sponse properties of one-dimensional systems. Here we 
show that exciton confinement dictates the evolution 
(saturation) of the optical response of one-dimensional 
chains as a function of the chain length. Therefore the 
failure of EXX, available local or semilocal or CDFT 
approximations is related to their inability to describe 
strong excitonic effects in anisotropic low-dimensional 
systems. We show that the xc-kernel has a hyper-non lo- 
cal behavior and memory dependence that is at least one 
order of magnitude stronger than in solids. The static 
and dynamical polarizabilities are both well described in 
simple H2 or in more sophisticated trans-polyacetylene 
molecular chains. 

The polarizability of low-dimensional systems is pro- 
portional to the spatial average of the reducible polariza- 
tion function, x (r, r';ui) that can be obtained from the 
solution of the TDDFT equation 



x M = xo M 1 + Ihxc M x M 



The exchange-correlation part of the kernel Jhxc 



(1) 



2 



f Hartree 

+ f^n as derived in Refs. [H, [lj] mimics to 
a good extent the MBPT results for the polarization 
function based on the solution of the Bethe-Salpetcr (BS) 
equation for the electron-hole dynamics in the basis of 
the Go Wo self-energy quasiparticle states [HI- This 
kernel is first order in the screened coulomb potential 
fxc 1 = Xo 1 Wxq 1 , where xo is the non-interacting 
quasi-particle response function, includes the cor- 



rect behavior lim q ^ fj^ K (y,q ) ~ — (7 + flu 2 ) /|q| 2 in 



the long-range regime The two constants, 

7 and P, measure respectively the degree of spatial non 
locality and memory effects. In spite of its simplicity and 
in contrast to EXX and CDFT, this fully captures 
the excitonic features in the absorption spectra of bulk 
semiconductors, insulators as well as surfaces Q. 
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FIG. 1: (Color on-line) Axial optical polarizability spectra of the 
molecular H2 chain A, (Peierls distorted, left frame) and chain B 
(nearly equally spaced H-atoms, right frame). While the present 
fxc kernel yields an excellent agreement with the Many-Body cal- 
culations in both cases, ALDA works only in the small gap case 
when the exciton is not localized. The excitonic wavefunctions 1 1 SH 
corresponding to the main absorption peak are shown fixing the the 
hole in the middle of the central hydrogen molecular bond. The 
delocalisation of the exciton wavefunction for the very small gap 
(chain B) is connected with a successful description of the spectra 
by the ALDA (see text). 

In order to illustrate how the kernel f^f- solves the 
failure of ALDA (and of any semi-local approximations) 
we consider two Peierls distorted infinite molecular H2 
chains . In chain A (see Fig. [1]) the inter-molecular 
distance is set to 2.5 a.u. and the intra-molecular dis- 
tance is 2.0 a.u. (as in previous works 0, S, El). This 
chain is a semiconductor with an LDA gap of 2.28 eV. 



In chain B the inter-molecular distance is reduced to 
2.05 a.u. and the LDA gap gets very small (0.26 eV). 
In Fig. [T] we compare the calculated TDDFT absorption 
spectra for chains A and B in the ALDA and MBPT 
f^. As the electronic density is more homogeneous in 
chain B than in chain A the ALDA gives an almost indis- 
tinguishable dynamical polarizability if compared to the 
BS calculations (see Fig [T] right panel). For the chain A, 
however, the ALDA performance worsens considerably, 
and the main absorption peak is not at all reproduced. 
The kernel, instead, yields a dynamical polarizabil- 
ity in both chains almost indistinguishable from the BS 
calculations. Similarly to the case of wide gap insula- 
tors [l6l | the correct description of excitonic states require 
the j^c" kernel to have non local Fourier components. 
This means that (tu, q) is a matrix, whose size (a 
few hundred reciprocal space vectors in the present case) 
is intimately related to the localisation of the excitonic 
state 0. 

The failure of the ALDA can be understood by look- 
ing at the spatial dispersion of the excitonic state cor- 
responding to the main absorption peaks shown at the 
bottom of Fig. [1] for both chains. In chain A, where 
ALDA does not work, the exciton is confined within few 
i?2-nionomers (~ 36 a.u.). Consequently, the excitonic 
dispersion introduces a characteristic length given by the 
exciton linear extension. In chain B the exciton is ba- 
sically spread all over the chain. If we now look at the 
xc-kernel for the two chains we see very drastic dif- 
ferences: whereas in chain B both 7 and (3 are almost 
zero, in chain A 7 - 16.92 and f3 ~ 0.45 eV 2 [23j], which 
is more than one order of magnitude larger than in solids, 
where 7 - 0.1 - 1.5, while (3 ~ 1 - 34 x 10" 3 eV" 2 [B]. 
This large value of (3 reflects the dynamical rcnormaliza- 
tion of the excitonic energy, due to memory effects. More 
importantly the static limit of the long-range part of the 
total TDDFT kernel, linwo f Hxc (w, q) ~ (4tt - 7 )/|q| 2 
is negative. This means that the spatial non locality of 
is stronger than the Hartree term (24j. This hyper 
non locality and the memory effects in cannot be 
captured by any static, local or semi-local approxima- 
tion. 

Turning now to finite-size effects, we show in Fig. [2] 
the results for the static polarizability of finite length H2 
chains within different approximations for the xc-kernel. 
Both ALDA and CDFT yields a very slow saturation of 
a as a function of the chain length and a strong overesti- 
mation of the static polarizability compared to accurate 
quantum chemical coupled cluster results CCSD(T). The 
TDDFT results obtained using the kernel partially 
reduce the ALDA and CDFT overestimation, showing a 
faster saturation. However the agreement with CCSD(T) 
is not yet satisfactory. We have traced the residual 
discrepancy to the lack of SIC in the LDA wavefunc- 
tions used to construct the xc-kernel. This can be eas- 
ily corrected by recomputing the kernel, using 
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Hartree-Fock (HF) self-intcraction-free wavefunctions to 
build xo- 
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FIG. 2: (color on-line) Axial polarizability per monomer of the 
molecular H2 chains. TDDFT calculations within the ALDA, 
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LR+SIC 



(see text) are compared with coupled clus- 



ter CCSD(T) results 0. and CDFT using the VK kernel [lj. The 
real-space excitonic wavefunctions of selected chains, with the hole 
placed in the middle of the central hydrogen molecule bond, are 
also showed. 

This new kernel yields axial polarizabilities in excellent 
agreement with the CCSD(T) results. The SIC increases 
the xc-kerncl spatial non locality factor 7, which trans- 
lates into further confinement of the excitonic states, 
compared to the LDA-bascd results. This can be 
rationalized by looking at the excitonic wave-functions 
shown in the bottom of Figj2] The probability of finding 
an electron at the end of the chain when the hole is placed 
in the middle of the bond of the central H2 molecule 
(that measures excitonic confinement) quickly decreases 
while increasing the chain length, being almost zero for 
for N = 9. The polarizability follows the same trend, sat- 
urating when the excitonic wavefunction dispersion does 
not change anymore. We should mention that, for small 
chains, ALDA gives good results as the excitonic wave- 
function is spread over the whole molecule, and the elec- 
tronic homogeneity is not strongly perturbed (like in the 



chain B case discussed above). It is clear that exciton 
confinement plays a major role in the determination of 
the response properties of Hi molecular chains providing 
microscopic support to the relevance of the ultranonlo- 
cality concept in exchange-correlation functionals. 
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FIG. 3: (color on-line) Calculated absorption spectra for the infi- 
nite (left frame) and polarizability per monomer for the finite (right 
frame) trans— polyacetylene chain. ALDA is compared with the 
fxc'- The quasi-particle gap and the position of the experimental 
absorption at 1.7 eV [27j are also showed. The axial polarizabil- 
ity for the finite chains calculated within HF |2SH and MP2 [j| are 
also reported. In the bottom we show the excitonic wave function 
corresponding to the bright 1.69 eV peak of the infinite chain. 

The H2 chain represents an extreme case of sys- 
tem with few electrons, where confinement is crucial, 
while electronic screening is negligible (this is the reason 
why time-dependent HF works fine). For these reason 
TDDFT kernels that work in the case of H2 chain may 
fail in more complex molecular systems where electronic 
screening becomes more important. This is the case of 
the trans-polyacetylene (PA) molecular chain. The PA 
chain has been extensively studied in the past within the 
BS scheme (25|. It has been showed that the main ab- 
sorption peak observed experimentally is excitonic, with 
0.44 eV binding energy. On the contrary all available 
calculations within TDDFT do not reproduce the exper- 
imental results: ALDA spectra (see Fig. |3|) is too low in 
energy and the main absorption peak shape is more simi- 
lar to a Van Hove singularity than to an isolated excitonic 
peak. Within an EXX calculation it is possible to obtain 
a good agreement with the experiment only by neglecting 
the x-kernel, i.e. in the independent particle approxima- 
which is clearly inconsistent. For finite length 



tion 

chains, as in polyethylene, EXX results are, again, very 
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similar to HF 0] that, in the case of PA largely over- 
estimates the quantum chemistry MP2 results (inset of 
Fig. [3]). As shown in Fig. O the present f£R kernel, in- 
stead, is in very good agreement with the MP2 results 
in the finite length case and with the experimental ab- 
sorption peak position in the infinite length case [26| . 
MP2 calculations [B| predicts the PA chain polarizability 
to saturate around f 5 monomers, in excellent agreement 
with the excitonic extension obtained with the present 
kernel (shown in the bottom of Fig. [2] for the bright 
PA exciton at 1.69 eV). As a result even in the more 
complex PA chain the polarizability saturation naturally 
correlates with the excitonic localisation length fully con- 
firming the physical picture that emerged from the pre- 
vious results for the H2 chains. 

In conclusion we have shown that TDDFT success- 
fully explain the optical saturation in molecular chains in 
terms of excitonic confinement. The Many-Body based 
xc-kernel, with the correct long-range tail and includ- 
ing, in the H2 chains, self-interaction corrections yields 
static and dynamical polarizabilities in excellent agree- 
ment with accurate quantum chemical calculations. We 
have proved that there exists a close link between the 
excitonic spatial extension and the axial polarizability in 
the TDDFT framework, giving a sound interpretation of 
the severe breakdown of the local-density approximation 
in anisotropic structures. 
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